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La estimación dela vulnerabilidad de os culivos agriolos 
al cambio climático se hace principalmente ignorando la 
¿iteración posible de la feilidad del suelo stibuibe al 
mismo cambio climático. E ojeivo principal del presente 
Arabaj fue estimar el papa dea altración dela feria del 
“suelo en ls predicciones ce La productividad del cultivo 
cd alo (Cua Arabic) ll el siglo XX para sl delas 
regonescafetaleras productoras as importantes del estado 
de Veracruz, México. Se han eunsiderado tres modelos 
de cculación global boo dos escenarios de foramiento 
aditivo. Se ha aplicado un modelo de desarllo del 
culivo en función de sus caacerísicas biológicas y las 
carateíicas climáticas propuesto. por JLASA/FAO y 
vtlizado ampliamente en el mundo. Los csfulos de la 
rodueividad deca cereza paca ls condiciones climáticas 
actuales y su comparación con los datos separados por 
la SAGARTA señalan que los rendimientos calculados 
son conflables. El coeficiente de comlición emir los 
"rendimientos callados y sbservados es igual a 093 con 
deror estándar de 48. A nal del siglo 20 e espera hasta 
un 34% de seducción enla productividad deca, vaiando 
Fentrnegionescafotaleras debido principalmente alcambial 
cambio en la prcipitaión y en menor grado, al incremento 
de La temperatura del sv La Ignorancia de la alteración 
de la fertilidad del sueo atsbuible al cambio climático 
posde causar erors e La estimación de la productividad 
de cl de hasta un 40%, por loque este facio no debe ser 
ecstimado. 


Palabras clave vulnerabilidad, agricultor, indice 
Fidroéraico local, indice integral de ferilidad del sueo, 
escenarios de cambio climático, estao de Veracruz, México. 


ISSN DIN74336« Treolaía y Cien del Agua, vo. V, má. 4 Juli agosto de 2 


Abstract 


Bráio, , Ni, Terms, L 6 Herr, $$. ue 
Agus, 20051. Exiate of he pct f Clinale Che on 
Sol Frey and Cafe Peach Vercrz, Mexico. Water 
Technology and Sence ln Sunis, 10.-116 


Estats o e ulnrally o op lo climate chungo pica 
Agron pose alteran dn Beer e ol ud by 
Als penso. The primary cc y e pee o us do 
stat e le alerts nl Ft e pss 
cof prluctin(of Aroca Lo en fe 205 contar 
A AS 
Mesta civcdatin mode ur uso bai 
fic cena, ea demon made rs pil srl 
o he Mig! and onto claras prposd y ASA 
TAO. Te mues extesiy used url he cala 
fo cf palacio used o cares ate condi nd 
Ae compis ni dta poi y SAGARPA indicada 
¡pe cl ae ea, coreatin con! Jete th 
callate src pis ás 135 tl standard ero 008 

By he el o 21 century, s much sa 30% elucti d cofi 
proc de sucia. hi tai aman of els prinurdy 
le canes nrtflmd Lose ente ur nar 
Hera. A lack of lag aut atera so fray 
cssal Ey ile change cn cate eros us his 4, e 
estimation o of prodactin. Tis factor hol hr ot de 
pe 


Keguords: Valora, sica local Seira nds. 
intogrl sol ferlty he, cate lago conri, tte 
VeneafMeso 


Recibido: 27/01/2014 
Aceptado; 06/03/2015 


5pp101 6 


ido eel, Estimación del impacto del cambio climático 


ye fertidad del sueo y prodocoidad 


Introducción 


Numerosos estudios recientes apuntan a que 
la variación del clima registrada en los últimos 
años es particularmente extraordinaria (Jacoby 
Se D'Arrigo, 1997; Mann, Bradley, de Hughes, 
1998, 1999; Rodney Se Freddy, 2002; Caballero, 
Lozano-García, Vázquez-Selem, E Ortega, 
2010). Se sabe de ciertos factores que pueden 
producir cambios en el clima, aunque no de 
manera precisa; mientras que algunos consi- 
eran que estamos ante un evento natural que 
+s parte de un ciclo, para la mayoría el cambio 
climático tiene su origen total o parcialmente 
en las actividades humanas, en particular en 
las emisiones de carbono a la atmósfera por 
el uso de combustibles fósiles y deforestación 
(Gonzalez et al, 2003; Magaña, 2004). 

Se considera que dicho cambio climático 
observado en la actualidad afectará en el futuro 
los recursos naturales de manera inevitable y, 
por la tanto, la productividad agrícola, inclu- 
endo el cultivo de café. Las investigaciones 
realizadas hasta el presente en el mundo y en 
México sobre vulnerabilidad de café al cambio 


climático consideran sólo el efecto directo de 
cambio de la radiación solar, temperatura dol 
aire y de precipitación, ignorando el efecto in- 
directo relacionado con la posible alteración de 
la frtlidad del suelo debido al mismo cambio 
climático, 

El estado de Veracruz ocupa el segundo 
Lugar nacional en producción de café, con una 
superficie de 153 mil hectáreas, que benefician 
de manera directa a wnas 300 mil familias de- 
dlicadas a este cultivo (Contreras, 2010), Hasta 
«el momento se han hecho varios estudios sobre 
vulnerabilidad del cultivo de café (principal- 
ente para Cofía rabicaL.) ante los escenarios 
de cambio climático esperado durante el siglo 
XXL La mayor parte de ellos representan análi- 
sis cualitativos (Camargo, 2010; Ville, Arizpe, 
Orellana, Conde, e Hernández, 2009; Haggar 
e Schepp, 2012), los cuales concluyen cómo se 
afectarán algunos procesos fisiológicos de desa 
srollo del cultivo, pero no permiten saber sobre 
el porcentaje del cambio en el rendimiento. 


Otros estudios estiman de manera cuantitati- 
va el impacto potencial del cambio climático 
sobre la productividad de café; estos estudios 
se basan en el uso de los modelos empíricos 
de tipo de regresión entre la productividad de 
café, en función de los factores climáticos y aún 
algunos factores económicos relacionados con 
el manejo de dicho cultivo (Fontagro, 2011; Ha- 
gar 8: Schepp, 2012; Laderach eb al, 2011; Lin, 
Perfecto, de Vandermeer, 2008; Paavola, 2008; 
Heakin, Gay, Estrada, e Conde, 2004; Jarami- 
lo ek al, 2011) o bien de los modelos teóricos 
(principalmente del ILASA-FAO, 2012) sobre 
el desarrollo del cultivo de café en función de 
sus características biológicas y las condiciones 
climáticas esperadas en transcurso del siglo XXI 
(Bunn eb al, 2013; Davis, Gole, Baena, de Moat, 
2012; Rivera:Silva e al, 2013) en los cuales han 
encontrado como resultado que se espera una 
pérdida en la cosecha de café (Coféa arabica L) 
en el estado de Veracruz en el transcurso del 
iglo XXI de orden de 10 a 25%, en función de 
los escenarios de cambio climático y modelos 
utilizados para tal análisis. 

No obstante, en dichos estudios se ignora el 
factor de la alteración potencial de la fertilidad 
del suelo atribuible al cambio climático a largo 
plazo. Durante los últimos años, se han hecho 
varios trabajos en México sobre la estimación 
cuantitativa de papel de cambio en la fertilidad 
del suelo y su efecto en el cambio dela producti- 
vidad agrícola (Nikolski, Bakhlacva, Contreras, 
¿e Ordaz, 2001; Contreras e a, 2002; Castillo et 
al., 2007; Terrazas, Nikolskii, Heerera, Castillo, 
de García, 2010); sin embargo, estos trabajos se 
dedicaron a estimar la vulnerabilidad al cambio. 
climático solamente de cultivos de trigo, maíz 
y frijol 

El objetivo del presente trabajo fue evaluar 
la vulnerabilidad de la producción del cultivo 
de café (Cofea ambica L.) en la zona centro 
del estado de Veracruz al final del siglo XXI, 
utilizando el modelo teórico de la HASA/ 
FAO (2012) y considerando no sólo el proceso 
fisiológico de desarrollo del cultivo en función 
de las condiciones climáticas; sino, además, el 
efecto del posible cambio en la fertilidad del 
suelo ocasionado porel cambio climático. 
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Materiales y métodos 


El presente estudio se llevó a cabo en la zona 
centro del estado de Veracruz, en 20 munici- 
pios que se encuentran distribuidos en las seis 
principales regiones cafotaleras del estado de 
Veracruz: Atzalan, Coatepec, Córdoba, Huatus- 
co, Misantla y Tezonapa. Dicha zona concentra 
90% de los productores y 93% del total de la 
superficie cafetalera de la entidad (Moguel de 
"Toledo, 1999; López, Díaz, $e Martínez, 2007) 
Las regiones se encuentran ubicadas entre 
latitudes de 18" 36' a 20* 4, y altitudes que van 
desde los 600 a 1.500 msnm; la precipitación 
medía anual es de 17173 mm; la temperatura 
media anual es de 202 “C, los tipos de suelo 
predominantes son andosoles, cambisoles, lito- 
soles y luvisolos, y las pendientes de los terrenos 
varían entre 0 y 30% (SMN, 2013; INIFAP, 2012), 
Las condiciones climáticas Kípicas para 
inicio del siglo XXI se obtuvieron utilizando 
las normales climatológicas para cada una de 
las estaciones climatológicas seleccionadas, 
mismas que se encuentran publicadas en la 
página web del Servicio Meteorológico Nacional 
(SMN, 2013). Las condiciones climáticas para 
final del siglo se obtuvieron de las 
nes desarrolladas por la Unidad Informática 
para las Ciencias Atmosféricas y Ambientales 
del Centro de Ciencias de la Atmósfera de la 


Universidad Nacional Autónoma de México 
(UNAM) hechas a partir de tres modelos de los 
15 modelos de circulación global disponibles en 
el proyecto de interoperación de modelos fase 5: 
MPLESMELR (Max-Plank Institute), GFDL.CM3 
(Geophysical Fluid Dynamies Laboratory) y 
HADGEM2-ES (Met Office Hadley) (Fernández, 
Zavala, te Romero, 2014), 

Se han considerado dos escenarios de 
forzamiento radíativo: RCP 45 (650 ppm de 
C02) y RCP 85 (1370 ppm); el forzamiento 
KCP 60 (720 ppm) no se ha analizado por 
considerarse un valor intermedio entre los 
seleccionados. 

Los tres modelos mencionados fueron utili 
zados tomando en cuenta las siguientes razones: 


— Aunque no existe gran diferencia entre los 
valores mensuales de temperatura y precip 
tación simulados por los diferentes modelos, 
los tres modelos utilizados se encuentran 
dentro del grupo de mejor resultado de 
desempeño, de acuerdo con la evaluación 
sealizada por Cavazos et al. (2013). 

— La resolución espacial de los 15 modelos 
presentados en el informe ina de la actua 
lización de escenarios de cambio climático 
para México, como parte de los productos 
dela Quinta Comunicación Nacional es de 
05" x05*(55 x55 km, aproximadamente) 
(Cavazos et al,, 2013). Sin embargo, los 
modelos utilizados fueron reestructurados 
mediante un proceso de reducción de 
vscala hasta 30" x 307 (926 x 926 m aproxi- 
madamento), tomando como referencia la 
climatología base 1950-2000 desarrollada 
por Hijmans, Cameron Para Jos Jrs (2005), 
quienes incorporan el efect topográfico y, 
porlo tant, se obtiene una mejor alternativa 
de distribución espacial delas variablos de 
cambio climático consideradas. Usar mode- 
los con estas características resulta de gran 
importancia para calcular rendimientos de 
café en regiones montañosas del estado de 
Veracruz, donde existe una gran variación 
climática entre zonas cafetales. 

— El propósito del presente trabajo fue estimar 
+l posible grado de alteración dela fertilidad 
de suelo y el cambio enla productividad de 
café, en función de los cambios esperados 
de temperatura y precipitación enel estado 
de Veracruz al final del siglo XXI. Para 
ello partimos del supuesto de que si los 
rendimientos calculados no van a cambiar 
die manera importante, esto significaría 
que puede ignorarse la vulnerabilidad de 
la producción de café al cambio climático; 
perosilos rendimientos dependen significa» 
tivamente de los cambios enla temperatura 
y precipitación simuladas por los modelos 
seleccionados, esto significaría que pos- 
teriormente será necesario profundizar y 
detallar este tipo de investigación. 
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Las características climáticas actuales de las 
regiones cafotaloras del estado de Veracruz se 
hicieron con estimaciones del impacto del cam- 
bio climático sobre la productividad de cereza 
de café (Cofiea Arabica L.) al final del siglo XXI y 
se presentan en el cuadro 1 

¡Como puede observarse, las precipitaciones 


varían desde una lámina de 9356 mm año" has- 
1a2755.97 mun; la temperatura media anual del 

Según Budyko (1973), Volobuev (197) 
Aydarov (1985), Nikolskii et. (2001), Contreras 
etal.(2002) Castillo ea. (2007), y Terrazas e 
(2010), las condiciones climáticas anuales para 
un año j pueden caracterizarse por el índice 
hidrotérmico local IHT: 


se encuentra en tn intervalo de 15 2 25"C, 


(0) 


er 


Donde Ru es la radiación neta promedio 
anual (MJ m* año"). es el calor latente de 
exaporación del agua (2512 M] mó men) y Pr 
esla precipitación media anual (mm año”), 

Los valores de JHT! son adimensionales 


y para México varían de 0-4 a 1 en las zonas 
tropicales húmedas; de 1 a 3 en las zonas 
semi-hámedas y semiáridas; y de 3 a 8 en las 
zonas áridas del norte del país. Como se ve 
en el cuadro 1, los valores del IAT%% para las 


segiones de referencia de producción de café.en 
el estado de Veracruz varían entre 0.73 y 195, 


lo que corresponde a sn clima semi-húmedo. 
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Estimación de la productividad agrícola del 
café (Coffe arabica L.) 


Los rendimientos del cultivo de café cereza se 
calcularon para los escenarios climáticos men- 
cionados anteriormente y para las condiciones 
«climáticas de inicio de este sigl, utilizando una 
ecuación propuesta por la HASA/FAO (2012) 
Yly= YL INSH! E] o 

Donde Y,,, es la productividad agrícola del 
cultivo de café (CofaarabicaL. ten kg ha" año" 
de café cereza) en los sitios de referencia para el 
escenario base inicio del siglo XXI ( = 2000) y 
para elfinal del sigo XXI (20752099) Y, es 
+l rendimiento potencial máximo o agraclimático 
de la materia seca económicamente aprovecha- 
ble que pueden producir plantas sanas con un 
suministro adecuado de agua y nutrimentos (en. 
kg ha? año"), para el escenario base a inicios 
del siglo XXI (¡=2000) y para el final de siglo 
XXI (¡ = 2075-2099); E) = índice integral de la 
fertilidad del suelo para el escenario base a ini- 
cios del siglo XXI = 2000) y para para el final 
del siglo XXI) = 2075-2099), el cual se supone 
puede ser diferente al final del siglo debido a 
ierta dependencia de la fertilidad del suelo en 
función de las condiciones climáticas (Contreras 
tal 2002; Castillo e? al, 2007; Terrazas et al, 
2010). 

Desde luego, el modelo no considera un 
efecto potencial de plagas o enfermedades 
(las cuales se supone también dependen del 
cambio climático) debido a que los modelos 
matemáticos existentes utilizados en las investi- 
gaciones para pronosticar elimpacto del cambio 
climático sobre enfermedades y/o plagas son 
empíricos (Van der Vossen, 2005). Es decir, se 
han desarrollado con base en observaciones 
hechas en algunos sitios de referencia, por lo 
cual no son representativos para su aplicación 
en otras regiones y zonas climáticas; además, 
tales modelos contienen una gran cantidad de 
parámetros empíricos, que no se sabría cómo 
aplicar para el caso del estado Veracruz. Por 


tal razón, dichos modelos no se utilizan en la 
presente investigación. 


Estimación de la productividad potencial 
UN) 

El cálculo dela biomasa y el rendimiento poten- 
cial o máximo Y/,, (en kg ha” año") del cultivo 
se realizó utilizando un modelo basado en prin- 
cipios ecofisiológicos (HASA/FAO, 2012). 


Yo =Bnric 5) 


Donde Bn es la biomasa neta de la materia 
seca total (en kg ha") e ICos el índice de cosecha, 
la fracción de Bn correspondiente al producto 
agrícola (dimensional). El valor de Bh (en kg 
ha“) se calcula con la siguiente ecuación: 


036b,L m5] 


Maat 


Dondeb,, esla tasa máxima de producción 
de biomasa bruta de un cultivo de referencia 
con el Índice de Área Foliar (IAF) igual a 5 (en 
kg ha! día") b,, depende principalmente de 
la radiación fotosintéticamente activa y de la 
concentración de CO, en la atmósfera; Les la 
fracción de la tasa máxima de crecimiento del 
cultivo, por cobertura incompleta del terreno 
cuando el índice de área foliar 1AF es menor de 
5: IAF <S (adimensional) 1 esla duración del 
ciclo normal del cultivo (en días) para el cultivo 
die café (Cofa ArabicaL., la duración de ciclo 
es de 270 días (LASA/FAO, 2012); C, esla frac» 
ción de la tasa de producción de biomasa bruta 
(como CH,O) que se pierde por la respiración 
de mantenimiento, misma que depende del tipo 
de cultivo (leguminosa o no leguminosa) yde la 
temperatura media del aire (kg día") 

Para la obtención de los valores de 5, se 
estimó la tasa máxima de producción de bio- 
masa, Pm (en kg CHO ha" Kv) mediante las 
expresiones siguientes: 


Para Pni'=20 kg CH,O ha" he: 


i 


dez015 pp. 
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2. =M(OS+0012,)p,+ (1) (0s+u0se Je (5) 


Para Pm <20 kg CH,O ha" hr 


»_=N(0820012,)t, +(1-N)(05+0.052, Je. 


(6) 


Donde N es un parámetro adimensional 
dependiente de la radiación fotosintéticamente 
activa (Ac) teórica o potencial diaria en el cielo 
totalmente despejado (en MJ m* día") y de la 
radiación global de onda corta (Rg, en MJ m* 
día'). Esta ecuación asume que la radiación 
fotosintóticamente activa (RFA) real esla mitad 
dela radiación global y que en un día nublado, 
la RFA es.el 20% de la teórica (Ac): 


Ac-05Rg 


() 


SAC 


b, es la tasa bruta de producción de materia 
seca para un cultivo de referencia hipotético (Kg 
CHO ha” ven días completamente nublados, 
«on dosel cubriendo por entero el terreno y una 
tasa máxima de producción de biomasa de 20 
kg ha*h Bes la tasa bruta de producción de 
materia seca para un cultivo de referencia hipo: 
tético (kg CHO ha” h) en días completamente 
despejados, con dosel cubriendo por entero el 
terreno y una tasa máxima de producción de 
biomasa de 20 kg ha". 

El factor de corrección por cobertura income 
plota del terreno (L) se calcula de la siguente 


SUAF<S 
L-O32A+0S05ILOgA(1AF) (8) 
SIAF25, entonces L=1 


Los valores de C, se pueden calcular con la 
siguiente ecuación: 


0) 


-C,(0.0044+0.00197.+090107" 


Donde T es la temperatura media mensual 
del aire durante el cilo del cultivo; C.jes la tasa 
de pérdida de producción de biomasa bruta 
por la respiración de mantenimiento a 30 “C 
(0.0108 para plantas no leguminosas, como en 
este caso). 

La información bibliográfica sobre 1C, N, JAF, 
ruta fotosintética, P., D,, h, y datos regionales 
sabre inicio y duración del ciclo del cultivo de 
café en México fueron obtenidos las publicacio- 
nes HLASA/FAO (2012) y De Wit (1965). 

En cuanto a la información para el cálculo 
según el periodo referido en los escenarios de 
cambio climático para el final del siglo XXI 
(=2075-2.099) se tomaron las siguientes con- 
sideraciones: El /C prácticamente no depende 
del cambio climático, por lo que no se consi- 
dera cambio; asimismo, el JAF se considera sin 
cambio (Cure de Acock, 1986; Anthony € Ziska, 
2000) la duración del ciclo del cultivo (N) se 
considera sin cambio, pues en la metodología 
o se consideran cambios para la fenología del 
cultivo de café debido a que no se cuenta con 


igaciones referentes al tema en nuestro. 
país; se asume que la transmisividad de la 
atmósfera no cambia de modo significativo 
debido a su propiedad diatérmica, por lo que 
los parámetros hr, y h permanecen sin cambio, 


Estimación del índice de disponibilidad de 
agua en el suelo (INSHI) 


El índice de Satisfacción de las Necesidades 
Hídricas se calculó a través del balance hídrico 
anual en la zona radical de tipo acumulativo; 
para elo se utilizó una escala mensual mediante 
la ecuación propuesta por Frtre y Popow (1986): 


ISNH] = ISNHL 


Donde ISNH) es el índice de satisfacción delas 
necesidades hídricas en el mes í en el año , es el 
ndice de Satisfacción de las Necesidades Hídricas 
para el mes anterior de ese mismo año; DEF, es 
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«el valor absoluto de la deficiencia de humedad 
en el mes ¡en el año ¿; Y NH! es la sumatoria 
dde las necesidades hídricas del cultivo en el año 
¡ desde el mes hasta el», considerando su cio 
(220 días, del 15 de fbrero a 15 de noviembre) en 
«el cual la planta tiene su producción de biomasa 
(LLASA/FAO, 2012) Para el cálculo delos valores 
del ISNH, y DEF, se ha considerado la evapo- 
transpiración del culivo (ETE) enel mos ¡del año, 
y la cantidad de agua quese infiltra en el suelo 
(nf)enclmes del año), 

Allen, Santos, Raes y Smith (2006) definen 
«como necesidad de agua del cultivo (o necesidad 
hídrico, NH) ala cantidad de agua requerida pa- 
ra compensar la pérdida por evapotranspiración 
del cultivo (ETE). Para el caso que nos ocupa, 
el cálculo de la evapotranspiración del cultivo 
(ETe) se realizó con base en sus recomendacio- 
nes, utilizando mínimos datos climatológicos; la 
evapotranspiración de referencia (ET) se calew- 
1ó por mes, utilizando el programa CROPWAT 
5.1(EAO, 2013) y posteriormente se multiplicó 
el valor obtenido por un valor cocficionte de 
transferencia Kc = 1.1 (Allen et al., 2006). 

La cantidad de agua infiltrada fue obtenida 
por diferencia entre la precipitación y el escurri- 
miento superficial por mes (SARH, 1991; CNA, 
2002); posteriormente se determinó la lámina 
de agua aprovechable almacenada y disponible 
para el cultivo con base en la textura de suelo y 
para una profundidad de 7D cm, a cual corres» 
ponde a la zona radical del cafeto. 

Desde lucgo existen cafetales que se encuen- 
tran bajo sombra, por lo que resulta necesario 
«considerar si existe un efecto de reducción de 
la evapotranspiración a través de los cálculos 
pertinentes con datos bibliográficos sobre distri 
bución de la radiación neta dentro o debajo dela 
apa delos árboles (Chang, 2001) Sin embargo, 
de acuerdo con observaciones experimentales 
sobr el cambio del microclima y de la produ 
tividad de café en el caso de la sombra (Siles, 
Harmand, de Vaast, 2010) se ha encontrado que 
si bien la temperatura de las hojas de cafeto se 
reduce, también existe una competencia por el 
agua del suelo, dado como resultado final una 
producción de café similar para las condiciones 


de cafetales a libre exposición y bajo sombra. 
Considerando que el objetivo principal de los 
cálculos era estimar el cambio relativo del rendi- 
miento al final del siglo XXI (122%) como una 
fracción del rendimiento al inicio de este siglo 
(012%), se estimaron los valores relativos ISNFP 
5200 [SNE Ep2S20m) Ea y pps pp 
a, concluyendo finalmente que se puede igno- 
tar el efecto de la sombra y asumiendo además 
que la arquitectura de la vegetación cambiará 
significativamente con el tiempo. 


Estimación del índice integral de fertilidad del 
suelo (E) 


Elúndice integral de fertilidad ageíola(F) para 
el escenario base se ha calculado según lafór- 
mula propuesta por Pegov y Jomyakow (1991) 
modificada. por Nikolski, Castillo-Álvarez, 
Bakhlacva, Román: Calleos y Maslow (2006) 


Donde Fes la fertilidad del suelo (o agrico- 
la) para el escenario base a inicios del siglo 01 
(¡=2000); MO, K y pH sonlos valores modales 
de contenido de materia orgánica en el suelo, 
contenido de fósforo y potasio disponibles 
para el cultivo y de pH típicos para los suelos 
dle cafetales con el mismo índice climático IHT 


para el escenario base a inicios del siglo XXI 
(¡= 2.000); MO, P.,. y K,,,. son los valores 
máximos observados de MO, P, K en los terre- 
mos estudiados (/= 2.000) 

El valor E, es adimensional, varía entro 1 y 
0.1 y corresponde al suelo más fértil potencial» 
mente; O, al suelo completamente degradado e 
infértil, La ecuación (11) es útil porque utiliza la 
información limitada sobre algunas propiedades 
de suelos de la república mexicana disponible 
en los datos de INEGI (1988, 2004) e INIFAP 
012). 

La estimación de la alteración de la ferti- 


lidad de los suelos al cambio climático está 


i 


Brgido eel, Estimación el impacto del com 


lmátco obre frdad del sul y productividad de 


basada en la Ley Geográfica de Zonificación 
del Suelo (Budyko, 1974; Volobuyey, 1974) 
De acuerdo con esta Ley, los valores modales 
regionales de algunas propiedades químicas, 
físicas y biológicas de los suelos vírgenes, 
geomorfológicamente homogéneos (con la 
misma topografía, condiciones geológicas e 
hidrogeológicas, textura y mineralogía del 
subsuelo y tiempo de formación) dependen 
principalmente del índice climático IHIT típico 
para el escenario base a inicios del siglo XXI 
( = 2 000), y para el final del sigo XXI 
( - 2075 - 2 099) (Volobuyew, 1973; Aydarov, 
1985; Niko'skii el, 2006; Castillo e a, 2007), 

En el cálculo del 1H, para el cálculo de la 
radiación neta (R,), se empleó la ecuación (12) 
(Allen et al, 2006) 


R,=(I-a)R K, 10) 


Donde «es albedo o coeficiente de reflexión 
del cultivo (adimensional), con un valor de 0.15 
(Jaramillo £ Gómez, 1989; Jaramillo, 2005); R, y 
E,,son la radiación solar entrante (o radiación 
global) yla radiación neta de onda larga, respec» 
tivamente, en el año / (MJ m* día”) 
Después se realizó un interpolación espacial 
de los valores de JHT para el escenario base 
2.000) para toda la zona cafetalera utilizando 
el método de Kriging con el sofware AreGis 1.10 
(ESRI, 2010), 
Al establecer relación cuantitativa entre los 


utilizando las bases 


valores modales de E? 
existentes de los datos de la distribución actual 
delas propiedades de suelos vírgenes y de cafe= 
tales de México, en función del índice climático 
JHT2>%, se puede suponer que según la Ley 
Geográfica de Zonificación del Suelo, debería 
«conservarse dicha relación F29%(1217299) al final 
del siglo XXI, cuando los valores regional: 

del índice climático cambiarán hacia IHT295 
+, considerando que el cambio climático, en 
promedio anual, ocurre en forma bastante lenta 
como para conservar la relación F2W9(1/1720%) 
Entonces, al conocer los valores IFI7295=2%) para 
los sitios de referencia de producción de café, se 


puede estimar cambio del factor de fertilidad en 
estos sitios F20509, utilizando la relación en la 
gráfica FIOO(ATI0s) 

La explicación esquemática de la relación 
existente entre la fetlidad y el JET para presen- 
te y su posible cambio en un fututo se presenta 
enla figura L 

El procedimiento para establecer la relación 
12% (1EIT2%9) ha sido descrito en los trabajos 
de Nikolskii eta. (2001, 2006) y Contreras eta. 
(2002). Para cl presente trabajo incluye princi- 
palmente las etapas siguientes 


1. Selección y sistematización de sitos geográ- 
ficos y de los datos sobre las propiedades 
mencionadas de los suelos vírgenes y de 
cafetales de la república mexicana utilizan- 
do el conjunto delas cartas de INEGI (2004) 
e INIFAP (2012). En el presente trabajo, el 
procesamiento de la información edafoló- 
pica fue automatizado a través del soft 
ArcGis 10.0; para el análisis, se utilizó el 
continuo de elevaciones mexicano 3.0 (CEM 
3.0), cuya resolución es de 15 m, y utilizando 
las bases de datos georreferenciadas sobre 
las propiedades de suelos de los pozos eda- 
fológicos (INEGI, 2004, 2013; INIFAP, 2012) 

2. Procesamiento estadístico delos datos sobre 
las propiedades de los suelos y el índice de 
fertilidad, búsqueda de sus valores modales 
X y los intervalos de confianza X +20 (don- 
dle a corresponde a la desviación estándar 
de los logaritmos naturales de los valores 
parciales de cada propiedad) en cada uno 
de los intervalos establecidos en la escala 


del índice climático JHT*%, En las inves- 
tigaciones previas fue establecido que la 
distribución estadística de cada propiedad 
tiene carácter lognormal (Nikolskii e? al, 
2001; Contreras el al, 2002). 

3. Establecimiento de la(s) curva(s) F; 
(ET) de mejor ajuste a los valores 
modales de esta propiedad, utilizando el 
programa CurveExpert 1-4 (Hyams, 2010), 
marcando además los intervalos de confian- 


Figido eel, Estimación del mpaco del cambio ciáti 


ve fertilidad del suo y producida de cof en Vo 


m5 


mr 


E pa 


gar. Esquema de predicó del cambio del fic integral de frida del suelo en una zona cafeta de nerncia 
en fación del cambio del nic climático IT en tanscurso del siglo XI. La gráfica UT) corsporade aL dependencia 

smos valores mada regale de en funeión delos alos regionales promedio anuales de HT l nico delia XXL 

Ty ES coresponden a uns de referencia enel inicio de est siglo. Los pares INT ¿Sy EPI O IATI=O y pap 


corespenden alos scenaio ponibes de cambio del nic clmático (4T 209 o 1H. 
Final delsiglo XX y alos valores sspocivos del índice de fetidad de pueo (E 


+1) el miso sii demencia al 


oro, 


4.. Determinación del valor del índice climático. 
IHTEC0% en cada sitio de referencia y del 
valor del índice correspondiente de fertili- 
dedo 


No se considera el efecto directo del aumen- 
to de la concentración creciente de CO, en la 
“atmósfera a finales del siglo XXI sobre el cambio 
de contenido de materia orgánica y el índice de 
fertilidad debido a que no es significativo (Baz- 
zaz $e Sombrock, 199; Kn, Prentice, Holland, 
Ae House, 2005) 

Por último, la estimación del cambio relativo 
enel rendimiento del cultivo de caféal final del 
siglo XXI (AY, en Js), producido para los esce- 
arios seleccionados, se realizó con la siguiente 
ecuación: 


1J100 (13) 


"También se ha estimado el cambio relativo de 
los componentes del rendimiento en el modelo: 
AY... AINSH y AF, utilizado las siguientes 


(13) 


AY] 100 (15) 


i 


aL VI nm. 4 o 
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Resultados y discusión 


Con objeto de verificar la confiabilidad del 
modelo para pronosticar la productividad del 
café, se estableció una comparación entre la 
producción de café cereza promedio anual en el 
periodo 2003-2011 y media entre los municipios 
de estudio (SIAP, 2013), 

La ecuación de regresión obtenida con el 
programa Curve Expert 1.4 (Hyams, 2010) es: 


0074099 3 10) 


son los rendimientos 
observado y calculado, respectivamente, como 
valores promedio anuales al principio del 


siglo XXI en los stos de referencia del estado 
de Veracruz; Vil, =2 4711 kg ha! año" y 
Villa. 20595 kg ha! año 

La diferencia entre los valores Yi y 
YEVL,. puede explicarse por la ignorancia de 
varios factores en el proceso de cálculo de los 
rendimientos, como aspectos económicos y 
tecnológicos de manejo del cultivo, fuctuación 
diaria de las temperaturas del are, riego de 
plagas y enfermedades, ente otros, Sin embar- 
$, para nuestro análisis es importante compar 
zar los rendimientos relativos, como fracción 
de los rendimientos máximos al final e inicio 
del siglo XXI. Debido a que existe una buena 
correlación entre los valores relativos YZ2"/ 
YSL, y VIVE, (Gara 2), consideramos 
que las estimaciones sobre el cambio relativo 
di la productividad esperados en el futuro son 
confiables 

En la figura 3 se presenta la dependencia 
dt os valores modales del índice de fertilidad 
del suelo F%% en función del índice climáti- 


ME 


Figura 2. Comparación entr Ls facciones de rendimientos observados y stmados com el modelo ILASA/FAD para 
cuicos de clima del periodo 20052011 
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gara. Dependencia delos valonss modales del índice fetlidad del sueo en función del índice climático ica paca 
Jos mos de cafetales del estado de Veracruz y ls suelos víngenesabicados en diferents zonas climáticas de Ménco en 
tertenos planos y Ladera con Vegetación font y dientes rango dependientes dentro el insrvao de 0 30%, También se 
muestran lsIneralos de comilanza 3 Za dé los valores modales con 95 de confiabilidad. Los suelos de caes del estado de 
Veracruz tenen os índicselimáticos promedio males 03 JT» 1,9, 


spica para los suelos de cafetales 
del estado de Veracruz y los suelos vírgenes 
ubicados en diferentes zonas climáticas de 
México en terrenos planos y laderas con 
vegetación forestal y diferentes rangos de 
pendientes dentro del intervalo de 0 a 30%. 
Asimismo, se muestran los intervalos de con- 
fianza + 20 de los valores modales con 95 de 
confiabilidad. 

La comparación preliminar de forma de las 
gráficas F2%% (1HT*%%) en el ámbito regional 
para suelos de los cafetales del estado de 
Veracruz con las mismas gráficas para suelos 
vírgenes ubicados en terrenos planos y laderas 
con vegetación forestal y diferentes rangos 
de pendientes dentro del intervalo de O a 
30% mostró que estas gráficas son similares. 
La distribución de los valores del índice de 
fertilidad F2'%% para suelos de cafetales de 
Veracruz en función del IHTY está dentro de 


rangos de variabilidad de los valores de 
para suelos vírgenes. Por esta razón, enla figura 
3 se presenta una gráfica general de distribución. 
El objetivo de la presentación generalizada de 
la gráfica F2%% (1FIT3%) fue estimar con mayor. 
seguridad el cambio relativo del índice de 
fortilidad E,en suelos de cafetales de Veracruz, 
en función del cambio del índice climático IET 
alfinal del siglo XL 

Los suelos de cafetales estudiados del estado 
dle Veracruz tienen los índices climáticos prome- 
dio anuales 0.73 <IHT?"< 195, que de acuerdo 
con la gráfica mostrada (figura 3) corresponden 
alos suelos bastante fértiles, 

Con base en los escenarios de cambio 
climático, al final del siglo XXI se espera un 
incremento promedio anual en la temperatura 
de 11 a25% (de 2 4*C), una disminución de 
107% enla precipitación (de 100 a 250 mm año"), 
en función del escenario y el sitio de referencia, 


i 
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Cuadro 2 Cambio estimado del rendimiento del cultivo deca (Co aaica 1) 8.) y de ss componentes (4%... 


INS 


y 5) en ls sion cafeteras del entro del Etado Verscaz al final el lo J0 (on período de 2752099) considerando, 
escenarios de cambio climático desarnllado por los modelos GFDL.CM, HADGEM y MDI ESM LK parados 2 trayectoria de 


úensciento de a concentración de COD.en la atmósfera al ial de este slohasta 650 pan y 137) pp. 


ía ía 
Región | modelo 
casaca || osas | 290. ]a20sn] ar, as] ar, Jada 
COM |] 
EN EEE) 7 78 32 
Arsalán — | HabcEno CET 330 [590 | 1100] 2070 
MPLESML 0 ao 7 3 22007 
GDL 2 | uz 030 | 045 5565 | 2638 
Contepec:— [HADGENZ EE 19077070 2555 [3053 
MPLESM ER 20 | 20 anar | 099 [2403 | 287 
GEDLCMS 730 [706 [1996 820 5 | an [aros 3507 
Córdoba | mabcena | 19 | 05 | 0 [2027500 ] 590 [oa 10 2000] 072 
| 5 951 [2009/9307 750 082 70] 037 
cra | 0] 0 50] | 0 20 22 957207 
huso [rapcae | 77 07 032 0 09 0 33 30] 
mar] os | 17 77070] 9 | 6% | 109 [2] 03 
amas | 005] 200 0% ]35| 87 
isso | mancene | 185 a jen[zo | os | 5 | 000 055] 25 
CESTAS EA PE PEO EE II ET E EE 
aca | 93 | 1039 | 1655 | 3690 2005] 2040 | 21039 1605 2506 | 1631 
tezonapa - [mabcenz | 959 [029 05011090 | 3066 | 1056 | 1029 100 [207] 05 
seres | 78] 0.09 | 650 [1056] 60 | 12 [900 1] 


y un incremento en la radiación global de 4 
10% (de 200 2400 Mjm* año"), 

Enel cuadro2 se presentan los resultados es 
perados de cambio mediante la modelación con 
la información relativa al cambio climático; la 
información presenta las estimación de cambio 
relativo delos rendimientos (Y, )y de loscom- 
ponentes de la productividad (AY, , AINSH y 
AF.) del cultivo de café (Cofc arabica L) para 
finales del siglo XXI (periodo 2075-2099) 

El análisis de la información contenida enel 
cuadro 2 señala que el cambio estimado de la 
productividad de café (AY,,,) al final del siglo 
XXX depende del escenario de cambio climático 
y delas rayectoras de crecimiento dela concen 
tración de CO, en la atmósfera, Se espera que la 
productividad de cafí al final del siglo XXI sea 


menor ala actual, con pérdidas desde un 10% 
hasta 34%. Se puede notar que los resultados de 
estimación de vulnerabilidad de café al cambio 
climático pronosticado por diferentes modelos 
tienen diferencia por lo general entre 30%. 

La mayor vulnerabilidad del cultivo de café 
(Cofeaaraca 1) al cambio climático se espera 
en as regiones de Coatepec y Huatusco, donde 
se estima una disminución de la productividad 
AY de 25 a 30%, encontrándose que el mayor 
impacto sobre el cambio de la productividad 
corresponde a la deficiencia de agua para de- 
sarrollo del cultivo. El cambio esperado del ín- 
dice de satisfacción de las necesidades hi 
(AISNH) debido a disminución de preciit 
y crecimiento dela evapotranspiración en algu 
as regiones se estima de 10 15%, mientras que 
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la productividad máxima potencial relacionada 
con crecimiento de la temperatura y de la radis 
ción solar se espera menor que la actual en las 
mismas regiones de 2a 6%, 

El análisisseñala que el cambio climático va 
a causar una alteración de fertilidad del suelo; 
el índice integral de fertilidad (E) va a cambiar, 
situándose estos cambios entre disminución 
hasta 5% y aumento hasta 15%, en función de 
la región; concentración de CO, en la atmósfera 
y el clima correspondiente. 

El error potencial en predicción del cambio 
de la productividad de café (gF) en caso de 
ignorar el factor de alteración de fertilidad del 
suelo debido al cambio climática se ha estimado 
con la fórmula siguientes 


td oo 15) 


Donde AYIÉ"**” (con F;) es el valor de cam- 
bio de la productividad estimada, tomando en. 
cuenta la fertilidad en el modelo y AY3%*" (sin 
E.) es el valor de cambio de la productividad 
estimada, sin tomar en cuenta la fertilidad en 


el modelo. 

Las estimaciones señalan que en caso de 
ignorar el factor de alteración de fertilidad del 
suelo E, debido al cambio climático se puede 
causar ervor en los cálculos de los rendimientos 
de café al final del siglo XXI hasta 50% o aun 
hasta cambio del signo de AY, lo que significa 
que este factor no debería ser ignorado. 


Conclusiones 


1. Se estima que la productividad de café 
cereza (Cofía orabica L) en la parte central 
del estado de Veracruz al final del siglo XXI 
va a reducirse entre un 10 y 34%, en función 
de la región cafotalera, trayectorias de cre- 
«cimiento de la concentración de CO, en la 
atmósfera y modelo que predice escenario 
del cambio climático, 


La mayor vulnerabilidad del cultivo de café 
al cambio climático se tiene en las regiones 
de Coatepec y Huatusco de ete estado. 

3. El análisis de vulnerabilidad de los come 
ponentes para la estimación de la produc» 
tividad del cultivo del caf: rendimiento 
potencial (Y) índice de satisfación de 
las necesidades hídricas (INSH), e índice 
integral de ferlidad (E), demostró que 
cada uno ejeres influencia en el valor total 
de la producción estimada y son impres- 
cindibles para determinar con un mayor 
grado de confiabilidad la variabilidad de 
la producción ante los escenarios de cambio 
climático. 

4. Laignorancia del factor de alteración de fe- 
úlidad del suelo debido al cambio climático 

puede causar ero en ls predicciones dela 

productividad de caf l final del siglo XXI 

de hasta-40%, lo que significa que ete factor 

no debe ser ignorado. 
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